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ABSTRACT 

We report X-ray studies of the environs of SN 1998bw and GRB 980425 using 
the Chandra X-Ray Observatory 1281 days after the GRB. Eight X-ray point 
sources were locahzed, three and five each in the original error boxes - SI and S2 
- assigned for variable X-ray counterparts to the GRB by BeppoSAX. The sum 
of the discrete X-ray sources plus continuous emission in S2 observed by CXO on 
day 1281 is within a factor of 1.5 of the maximum and the upper limits seen by 
BeppoSAX. We conclude that S2 is the sum of several variable sources that have 
not disappeared, and therefore is not associated with the GRB. Within SI, clear 
evidence is seen for a decline of approximately a factor of 12 between day 200 and 
day 1281. One of the sources in SI, Sla, is coincident with the well-determined 
radio location of SN 1998bw, and is certainly the remnant of that explosion. The 
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nature of the other sources is also discussed. Combining our observation of the su- 
pernova with others of the GRB afterglow, a smooth X-ray light curve, spanning 
~ 1300 days, is obtained by assuming the burst and supernova were coincident at 
35.6 Mpc. When this X-ray light curve is compared with those of the X-ray "af- 
terglows" of ordinary GRBs, X-ray Flashes, and ordinary supernovae, evidence 
emerges for at least two classes of lightcurves, perhaps bounding a continuum. 
By three to ten years, all these phenomena seem to converge on a common X-ray 
luminosity, possibly indicative of the supernova underlying them all. This con- 
vergence strengthens the conclusion that SN 1998bw and GRB 980425 took place 
in the same object. One possible explanation for the two classes is a (nearly) 
standard GRB observed at different angles, in which case X-ray afterglows with 
intermediate luminosities should eventually be discovered. Finally, we comment 
on the contribution of GRBs to the ULX source population. 

Subject headings: gamma-ray bursts 

1. Introduction 

One of the most exciting developments in the study of gamma-ray bursts (GRBs) was 
the discovery, in 1998, of a GRB apparently in coincidence with a very unusual supernova 
of Type Ic (Galama et al. 1998). This coincidence of SN 1998bw and GRB 980425 offered 
compeUing evidence that GRBs are indeed associated with the deaths of massive stars, 
and that, at least in some cases, GRBs go hand in hand with stellar explosions (Woosley 
1993; MacFadyen & Woosley 1999). The large energy release inferred for the supernova also 
suggested a novel class of explosions, called by some "hypernovae" (Paczynski 1998), having 
unusual properties of energy, asymmetry, and relativistic ejecta. 

However, this identification was challenged on two grounds. First, there were two vari- 
able X-ray sources identified with BeppoSAX in the initial 8.0' radius GRB error box; one 
was not the supernova. Second, if it were associated with the nearby supernova, GRB 980425 
was a most unusual burst with gamma-ray energy per solid angle roughly four orders of mag- 
nitude less than typical. Further, the BeppoSAX decay of the SN 1998bw associated X-ray 
source was much slower than the typical GRB X-ray afterglow decay (Pian et al. 2000). 
Interestingly, observations (Pian et al. 2003) of the two sources in March 2002 with the 
X-ray Multi- Mirror (XMM) telescope revealed that the X-ray emission of the SN 1998bw 
associated BeppoSAX error box had decreased at a faster pace than expected by a simple 
extrapolation of the earlier measurements. Moreover, the second XMM source was found 
to consist of a number of faint point-like sources, whose integrated emission was consistent 
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with the average brightness measured with BeppoSAX (Pian et al. 2000). 

Strong support for the GRB-SN association came from the very recent spectroscopic 
detection (Hjorth et al. 2003; Stanek et al. 2003) of a supernova (SN 2003dh) in the optical 
hght curve of GRB 030329. In this case, the SN detection was obscured by the extreme 
optical brightness of the GRB afterglow. While the models of MacFadyen and Woosley 
(1999) suggest that there can be adequate ^^Ni to make the GRB-supernova bright, it is 
now becoming evident that we must also account for the potentially much brighter optical 
afterglow of the GRB itself. 

Here we report the results of a study of the environs of SN 1998bw and GRB 980425 
using the Chandra X-Ray Observatory (CXO) 1281 days after the GRB. This study had 
several goals. First, given the intervening three years, has the evidence strengthened for the 
GRB-SN association? We believe that it has (§ 4.1). Second, how does the X-ray light curve 
of GRB 980425, measured across 1300 days, compare with those of other GRBs and with 
other kinds of high energy transients - in particular X-ray flashes (XRFs; Heise 2003) and 
supernovae? What does the comparison tell us about the nature and origin of GRBs? We 
find that it provides evidence for a common theme underlying all these events - a powerful 
asymmetric supernova with relativistic ejecta along its polar axes and an observable event 
that varies depending upon the viewer's polar angle (§ 4 and § 5). 

Finally, we are interested in the environs of SN 1998bw. Aside from the one supernova, 
does this region show evidence for unusual stellar activity as might characterize a vigorously 
star forming region (§ 3.1)? SN 1998bw offers the best opportunity to study a GRB site up 
close, and one should take every advantage of that. 

2. Prompt Observations of GRB 980425 and SN 1998bw 

GRB 980425 triggered the Burst And Transient Source Experiment {BATSE ) on board 
NASA's Compton Gamma-Ray Observatory (CGRO) on 1998 April 25, 21:49:09 UT; the 
event was simultaneously detected by the BeppoSAX Gennma-Teiy Burst Monitor (GRBM) 
and Wide Field Camera ( WFC) . The burst consisted of a single peak of ~ 23 s duration, with 
peak flux and fiuence (24-1820 keV) of (3.0 ± 0.3) x 10"^ erg/cm^ s, and (4.4 ± 0.4) x 10"^ 
erg/cm^, respectively. Galama et al. (1998) observed the WFC 8.0' error box with the New 
Technology Telescope (NTT) at the European Southern Observatory (ESO) on April 28.4 
and May 1.3 UT and, in the error box of GRB 980425, they found supernova SN 1998bw, 
located in an HIT region in a spiral arm of the face-on barred spiral galaxy ESO 184— G82, 
at z=0.0085, corresponding to a distance of 38.5 Mpc (Galama et al. assumed a Hubble 
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constant, Hq — 65 km/s Mpc, but in the following we use Hq — 72 ± 8.0 km/s Mpc as 
measured by Preedman et al. (2003), placing SN 1998bw at 35.6 Mpc). 

On 1998 April 26-28, ten hours after the GRB, the BeppoSAX Narrow Field Instruments 
(NFI) observed (Pian et al. 2000) the WFC error box and revealed two previously unknown, 
weak X-ray sources, ISAX J1935.0-5248 (hereafter SI) and ISAX J1935.3-5252 (hereafter 
S2) with an uncertainty radius of 1.5' each. SI included the SN 1998bw, but S2 was 4.5' 
away. Both sources were observed two more times with the A^F/ resulting, in the case of S2, 
in two detections and two upper limits. In contrast, during the six month interval spanned 
by all NFI observations, the flux, F, of SI followed a power-law temporal decay (Pian et 
al. 2000): F2-iokeV = 4.3(±0.5) x 10-^2(t/lday)"°-2 erg/cm^ s, a much flatter trend than the 
one observed for other GRB X-ray afterglows, but a decaying trend nevertheless. 

At 35.6 Mpc the apparent isotropic energy of GRB 980425 (7 x 10*^'' erg) was about four 
orders of magnitude smaller than that of 'normal' GRBs (Bloom et al. 2003). Moreover, 
independent of its connection with a GRB, SN 1998bw was extraordinary in many ways. Its 
light curve resembled a Type la supernova in brilliance, but the spectrum was more like Type 
Ic (H, He, and Si lines absent, but the spectrum was peculiar even for Ic). Thus, the two 
phenomena together presented a very interesting scientific puzzle whose solution required 
the combined superb resolution of the Hubble Space Telescope (HST) and CXO. 

3. CXO Observations at Day 1281 

CXO observed SI and S2 on 2001 October 27, for a total time on source of 47.7 ks. SI 
fell completely on ACIS-S3 (a back-illuminated CCD) and S2 only partially, with most of 
the error region falling on ACIS-S2 (a front-illuminated CCD), both operating in Time Ex- 
posure (TE) mode. The data were processed using the CIAO (v3.0.1)^ and CALDB(v2.34)2 
software. More specifically, we used the CIAO tool acis-process- events to ensure that the 
latest gain corrections were applied (those corresponding to our observation date). Further, 
we removed the standard pixel randomization, applied CTl corrections, filtered the data to 
include events with ASCA grades = 0,2,3,4, and 6, and applied standard GTls. We cor- 
rected the systematic offset in the aspect using the fix- ojf set thread.^; this results in an offset 
of 5RA = 0.16" and 6DEC = 0.43". To improve the ACIS-S3 spatial resolution we use 



^ http :/ /asc . harvard . edu/ciao/ 
^http://asc. harvard.edu/caldb/ 
^http://cxc.harvard.edu/cal/ASPECT 
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the method described by Mori et al. (2001) and Tsunemi et al. (2001) to adjust the event 
locations. These studies conclude that by defining the event location on a pixel as a function 
of the event grade (rather than placing the event at the pixel center) one can achieve ~ 10% 
improvements in spatial resolution. During our observation the X-ray background increased 
by ~ 50% with variations on time scales of a few kiloseconds. Since the point source emission 
is not significantly masked by such slow background variations, we chose to include all the 
data in our analysis to preserve our (limited) source counts. 

3.1. SI and S2: Source Identification, Locations and Energetics 

We used the source-finding method originally described in Swartz et al. (2002), accepting 
as detections all sources with a minimum signal-to-noise ratio (S/N) of 2.6. For source 
detection purposes wc searched images consisting of data between 0.3 — 8.0 keV to avoid the 
ACTS high energy background. We discuss below the sources within the 1.5' radius (la) NFI 
error circles of SI and S2 only. 

3.1.1. S2 

S2 was resolved into five sources (S2abcde; Figure Ic and Table 1). Given the very 
limited count number per source, we combined four of them in one spectral fit with a power 
law function and a hydrogen column density, iVn = 3.95 x 10^° cm~^ (corresponing to the 
Galactic absorption in the line of sight to GRB 980425; Schlegel, Finkbeiner & Davis 1998). 
Here we assume a spectral similarity between these sources to mitigate the difficulty of 
fitting individual spectra of very few counts each. This order of magnitude approximation is 
acceptable within our source statistics (see also Figure 2, where rudimentary count spectra 
are presented for Slabc, and S2cc). However, one source (S2d), fell into the gap between 
the ACIS CCDs S2 and S3 and its total counts had to be adjusted taking into account 
the reduced gap exposure time. Using the spectral index of -1 from the fit, we applied a 
conversion factor of 1.89 x 10"^-*^ erg/cm^ s (0.3 — 10.0 keV) per count/sec (0.3 — 8.0 keV) 
to the S2d counts and finally estimated the total (combined) flux in S2 (0.3 — 10 keV) to be 
3.0(3) X 10~^^ erg/cm^ s (here and in Table 1 the numbers in parentheses correspond to the 
1(7 errors in the last digit). 

The BeppoSAX flux for S2 was, however, calculated for a much larger extraction radius 
(3') to account for the extreme faintness of the detection within a 1.5' radius, which was at 
the level of the A^F/ confusion limit (Plan et al. 2000). Only ~ 60% of this larger error circle 
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is covered with CXO. Using the same source-finding algorithim criterion described above, we 
found a total of 18 sources within the enlarged area, with a total flux of S2{3') = 1.6 x 10~^^ 
erg cm~^ s. This is roughly the flux to be compared to the BeppoSAX S2 value at day 1. 
Subsequently we re-calculated the total BeppoSAX Hux within S2 at day 1, using the same 
spectral function derived with CXO and found it to be 4 x 10^^^ erg/cm^ s (with < 2a 
significance) . Taking into account the partial coverage with CXO of the NFI error circle of 
S2, we estimate that the day 1281 CXO fiux value is within a factor of 1.5 of the 3a detection 
limit of BeppoSAX, indicating that at best there was no siginificant variation of the sources 
within S2 over the last 3.5 years. This result is consistent with the XMM observations of S2 
(Pian et al. 2003). We further discuss the evolution of the light curve of S2 in § 4.1. 



3.1.2. SI 

We initially identified two sources 36" apart within the SI (1.5') error region"', one of 
which coincided with the location of SN 1998bw. However, further inspection of the latter 
resolved this source into two with a radial separation of ~ 1.5" (Figure lb). We fitted these 
sources simultaneously with two 2-D circular Gaussians to better estimate their centroids 
and found that they are both consistent with point sources. Hereafter we designate the 
X-ray sources detected within the SI error region as Sla (corresponding to SN 1998bw), 
Sib, and Sic. We have fitted a power law to each unbinncd non-background subtracted 
source spectrum, assuming the same A^h as for S2. Our fit parameters are derived using the 
C-statistic (Cash 1979), appropriate for low count data (Table 1). 

Further, we reanalyzed all archival HST and ESO/Very Large Telescope (VLT) obser- 
vations of ESO 184— G82, the host galaxy of SN 1998bw, concentrating on the immediate 
environment of the supernova. The small field of view of HST/ STIS (~ 50 ") contains only 
the sources Sla and Sib. In order to obtain accurate astrometry we, therefore, registered 
both the CXO and iJST images independently to an R-band observation obtained at the VLT 
on 1999 April 18. The 3.4' x 3.4' field of the VLT image contained three additional X-ray 
sources with apparent optical counterparts which we used to align the two fields. Finally, 
we aligned the HST and VLT images using eight, non-saturated point sources present in 
each image (both alignments were performed using IRAF and the tasks GEOMAP & GEOXY- 
tran). We were then able to project the relative position of Sla and Sib onto the i^^T field 



■^Here we are only considering the area corresponding to the NFI half power radius of 1.5', as most of the 
BeppoSAX signal for SI at day 1 was within this area (Pian et al. 2000). An enlarged radius region (3') for 
SI results in a total flux of 2.7 x 10~^^ erg cm~^ s, almost half the BeppoSAX value at day 1 
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with a positional accuracy of ~ This is shown schematically in Figures le,f, where we 
zoom in the region around Sla and Sib for both the CXO and HST data. It is obvious in 
the HST image that while there is a clearly identified optical counterpart to the SN 1998bw, 
there is no variable counterpart within the ~ 0'.'3 CXO error circle (la radius) of Sib; we do 
however, see a variable source just outside the error region. Levan et al. (in preparation) 
present a detailed study of this source together with three more transients within a radius of 
6" of the SN 1998bw as well as narrow field spectroscopy of the SN 1998bw environment. A 
counterpart search for the other six sources in SI and S2 in all available catalog data^ failed 
to identify any known objects at their positions. 



Table 1: CXO Sources Within SI and S2 



ID 


RA 


DEC 


Counts" 


S/N 


Index^ 


Flux<= 


Luminosity'^''^ 




(hms) 


(° ' ") 








erg cm~^ s~^ 


erg s-i 


Sla 


19 35 3.31 


-52 50 44.8 


24 


4.53 


1.0(3) 


8(2) X 10-^5 


1.2 X 10^9 


Sib 


19 35 3.23 


-52 50 43.4 


52 


6.47 


1.6(2) 


9(1) X 10"^^ 


1.4 X 10^9 


Sic 


19 35 0.56 


-52 50 21.2 


27 


4.69 


0.5(3) 


1.2(2) X 10"^^ 




S2a^ 


19 35 25.75 


-52 54 19.0 


10 


2.75 




4(1) X 10-15 




S2b^ 


19 35 25.22 


-52 54 21.2 


7 


2.66 




3(1) X 10-15 




S2c 


19 35 24.65 


-52 54 40.7 


22 


4.46 


1.7(4) 


5(1) X 10-15 




S2d-^ 


19 35 23.15 


-52 53 5.6 


9 


2.82 




8(3) X 10-15 




S2e 


19 35 14.03 


-52 53 51.9 


23 


4.19 


0.7(3) 


1.0(2) X 10-1^ 





"between 0.3 — 8.0 keV, ''for an energy spectrum of cx , "^between 0.3 — 10.0 keV, ''assuming the source 
is in ESO 184— G(S^ at 35.6 Mpc, '^fluxes for sources with less than 20 counts are derived using a PL index 
of -1.0, and a conversion factor of 1.89 x 10^^^ ergs s~^ cm~^ (0.3 — 10.0 keV) per 1 count/s (0.3 — 8.0 keV), 
■''source flux corrected for decreased exposure 

From Table 1, wc notice that the X-ray luminosities of Sla and Sib marginally exceed 
the Ultra Luminous X-ray source (ULX) threshold luminosity of L = 1 x 10^^ erg/s (Fabbiano 
1989).^ Further, inspection of the Sib light curve in Figure 2 reveals that the source was 
'on' during the first 20 ks of our CXO observation, remained 'off' for the following ~ 28 ks 
and possibly turned on again during the last 2 ks of the observation. Counting the GRB-SN 
source also as a ULX, we then have two of these sources in ESO 1M—G82. What is the 



^http: / /heasarc. gsfc.nasa.gov/ db-perl/W3Browse/w3browse.pl 

^Here we are not discussing the nature of Sic as there is no clear evidence for the association of this 
source with ESO 184-G82 
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probability that Sib resides in ESO 184— how common are ULXs in galaxies and what 
is the nature of Sib, this ultraluminous X-ray transient? 

To address the first question we use the log A'"-log S relation for the hard X-ray sources 
by Giacconi et al. (2001): A^(> S) = 1200 x ( ^xw-^^ )-i±o.2 gources/deg^ (the one for 
the soft sources is smaller, so we use this as a conservative upper limit). Substituting for 
5" = 9 X lO"-*^^ erg/cm^ s, we find A^(> S) = 400 sources/deg^. Correspondingly, we expect 
0.01 source within a radius of 10". We conclude that Sib is most likely associated with the 
spiral arm of ESO 184—082 and is not a background source. 

There are two extensive studies of ULX populations in normal galaxies (Roberts and 
Warwick 2000; Humphrey et al. 2003) based on data from Rosat and 0X0, respectively. 
While both studies find substantial evidence for a correlation between the ULX count number 
with star forming mass, the Rosat survey tends to underestimate their count numbers by 
a factor of 5 — 10 with respect to the 0X0. We have, therefore, used the Rosat luminosity 
distribution DN/dLss = (1.0 ± 0.2)LJ8^'^, normalized to a 5-band luminosity of lO^^L© to 
calculate the expected number of ULXs in ESO 184— G82, multiphed by a factor of 10 as 
indicated by the 0X0 survey. We derive a total corrected expected number of 0.28 ULXs in 
ESO 184—082, a factor of ~ 7 less than the actual observed number of 2. We conclude that 
if Sla and Sib are both in ESO 184—082, there is a somewhat unusual concentration of 
ULXs in this galaxy. Finally, ESO 184— 082 is one of 5 members of a galaxy group. In the 
0X0 data, which covers all members, we do not detect any other X-ray sources from this 
galaxy cluster. This practically means that, assuming that all members arc at the distance 
of ESO 184— 082, we detect only 2 ULXs from the whole cluster (one of which is a GRB/SN 
afterglow), both within a radius of 6". We discuss the results of narrow band photometry 
of ESO 184— together with the properties of other unusual transients found in the HST 
data in the SN 1998bw environment in detail in Levan et al. (in preparation). 

To test Sib for variability we extracted the source lightcurve in 1000 s bins and per- 
formed a Kolmogorov-Smirnov test against a constant flux distribution; this is the prefered 
method for the low number of counts (< 100) which are seen in the source. We derive 
Pks = 0.05 (where Pks is the probability that the source is constant); thus Sib is consistent 
with being variable at the 95% confidence level. The transient nature and the energetics 
of Sib indicate a possible microquasar similar to e.g., GRO J1915-I-105 in our own Galaxy 
(Mirabel & Rodriguez 1994; Greiner, Morgan & Remillard 1996). This superluminal source 
exhibits transient behavior in a wide variety of time scales and intensities, and has an ap- 
parent (isotropic) X-ray luminosity of ~ 7 x 10^^ erg s~^, well above its Eddington limit 
(Greiner, Cuby & McCaughrean 2001). Fabbiano et al. (2003) have found in a study with 
0X0 of the ULX sources in the Antennae galaxies that seven out of the nine ULXs are time 
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variable, most likely accreting compact X-ray binaries. Interestingly, the average co-added 
spectrum of the Antennae ULXs resembles that of Galactic microquasars (Zezas et al. 2002). 
We conclude that most likely Sib is a microquasar in ESO 184— 

4. Systematics of the X-ray light curves of High Energy Transients 

Figure 3 shows, on a common scale, all GRB afterglows with measurements covering 
several tens of days. There are, unfortunately, only a few curves because such measurements 
can only be made on GRBs that are relatively nearby, but their light curves should be 
illustrative. For GRB 030329, we used data points from both the Rossi X-Ray Timing 
Explorer (RXTE) and XMM reported by Tiengo et al. (2003). The original four points of 
GRB 980425 were reported by Pian eA, al. (2000); the fifth point is the CXO observation 
of this paper. For GRBs 021004, 010222, 000926, and 970228, we used data from Sako & 
Harrison (2002ab), Bjornsson (2002), Harrison aZ. (2001), and Costa et aZ. (1997), 

respectively. We reanalyzed all data, calculated the fluxes in the same energy interval (0.3—10 
keV) and then converted them into (equivalent isotropic) luminosities using a cosmology of 
VLm — 0.27, Qa = 0.73, and Hq = 72 km s~^ Mpc~^ to ehminate the distance dependence. As 
has been noted many times, GRB 980425 and its early afterglow falls orders of magnitude 
below the "ordinary" GRBs. 

We then compared these curves with both those of supernovae of all types and with 
XRFs. We collected, reanalyzed and converted flux data to luminosities as described above 
for the Type lib SN 1993 J (Uno et al. 2002; Swartz et al. 2002), the Type Ic supernovae 
19941 and 2002ap (Immler et al. 2002; Sutaria et al. 2003) and the Type II supernovae 
1998S and 1999cm (Poolcy et al. 2002). Each source is indicated with a different symbol in 
Figure 3. Finally we added the three XRF afterglow light curves available (011030, 020427, 
and 030723; Bloom et al. 2003; Butler et al. 2003) assuming a redshift ol z — 1 for each of 
them (there are no redshift measurements for any of these sources) . The slopes of the XRF 
afterglows agree well with those of the typical GRBs and their luminosities are comparable 
for the distances assumed. Recently, Soderberg et al. (2003b) have reported the counterpart 
identification of XRF 020903 at a redshift z = 0.251. To reflect this lower distance scale for 
XRFs, we plot on Figure 3 another set of light curves (dashed lines) corresponding to the 
luminosities all XRFs mentioned above would have, when placed at a distance of z = 0.251. 
They still fall well within the 'typical' GRB range; thus XRFs would have to be extremely 
nearby for their X-ray light curves to be distinct from the generic GRB X-ray afterglow. 

The resulting plot is striking in several ways. Despite the huge disparity in initial appear- 
ance, there are compelling indications of a common convergence of all classes of phenomena 
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plotted - GRBs, GRB 980425, XRFs, and the most energetic supernovae - to a common 
resting place, L ~ 10^^ — 10^° erg about three to ten years after the explosive event. 
GRB 980425, being the closest by far of any GRB ever studied, has the virtue of being fol- 
lowed all the way to the "burial ground" , but a simple logarithmic extrapolation of the GRB 
and XRF light curves places them squarely in this region as well. As the collapsar model 
has predicted (MacFadyen & Woosley 1999) and observations of SN 2003dh/GRB 030329 
have unambiguously confirmed for one case (Hjorth et al. 2003; Stanek et al. 2003), an 
energetic supernova is expected to underly all GRBs of the long soft variety (Kouvehotou et 
al. 1993). Zhang, Woosley & Heger (2003) have also predicted that a similar supernova will 
underly all XRFs. 

We proceed now to a wider GRB X-ray afterglow comparison with those plotted on 
Figure 3. We have often observed GRB afterglows with temporal decays which cannot be 
described by a single power law. Rather, these afterglows initially decay as or a bit 
steeper, and at later times the decay index becomes approximately —2. This steepening is 
attributed to the fact that the ultrarclativistic outflow is initially coUimated within some 
angle 9 (typically thought to be of order 10 degrees. Frail et al. 2001; see also Van Paradijs 
et al. 2000, and references therein). The transition to a steeper power law marks the time 
when the outflow begins to expand laterally, which occurs around the time when 9 ~ l/F, 
where F is the Lorentz factor of the blast wave. Typically, the 'break time' when this 
happens is around 0.3 — 3 days after the burst, but a few cases have much later breaks, if 
any (e.g., GRB 970508, 000418; Frail et al. 2001). If this is also true for the X-ray afterglows 
shown in Fig. 3, then they should be extrapolated to later times with a steeper power law 
than the average of the data shown. In that case, they would reach the luminosity level of 
SNe sooner, and after that evolve like the SNe, because the blast wave will have become 
non-relativistic. Thereafter, the system may simply evolve like an X-ray supernova with an 
energy that has been augmented with that of the initially relativistic blast wave. However, 
it is tantalizing that a few systems, e.g., GRB 030329 do not yet show such a steep decay 
even at 30-40 days, possibly indicating that other mechanisms than the standard coUimated 
afterglow contribute to the emission. For example, it has been proposed that the outflows 
of GRBs are structured (e.g., Meszaros et al. 1998, Rossi et al. 2002, Ramirez-Ruiz et al. 
2002), i.e., they eject material with ever lower Lorentz factors at ever larger angles from the 
jet axis, and this slower material affects the afterglow at ever later times, making it decay 
more slowly. In this context, it is interesting to note that in GRB 030329 there is evidence 
for two jet breaks, one at 0.5 days and one at about one week (Berger et al. 2003) from radio 
to X-ray data. Moreover, WSRT radio data at 1.4 — 5 GHz indicate that there is probably 
an even wider outflow of yet slower material (Rol et al. in preparation). It is not clear 
however, that this material produces enough X-ray emission to explain the slow late decays. 
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and so the cause of these may be altogether different. 

Prom Figure 3 alone, it is not clear whether ordinary GRBs and GRB 980425 are two 
distinct classes of events with different X-ray light curves, or they form the boundaries of 
a continuum of high energy transients that will eventually fill in the entire left side of the 
figure. A simple theoretical interpretation discussed in § 5 favors a continuum of events. As 
we discuss, at early times an off-axis observer sees a rising light curve, peaking when the jet 
Lorentz factor is ~ 1/^obsj and approaching that seen by an on-axis observer at late times. 
This leads us to believe that the low X-ray luminosity of the recently discovered GRB 031203 
may have also been due to its having being viewed substantially off-axis. 

4.1. The X-Ray Light Curve of Sources SI and S2 

Figure 4 summarizes the X-ray observations of sources SI and S2. The first four data 
points (or upper limits) up to day 200 are from BeppoSAX (Pian et al. 2000) and do not 
resolve individual sources within SI and S2; our CXO observations on day 1281 do resolve 
the sources. We consider two hypotheses: 1) that SN 1998bw and GRB 980425 were the 
same event, both happening within SI, and 2) that GRB 980425 occurred within the error 
box of S2 at a cosmological distance and was thus a more ordinary GRB. Hence the S2 
observations are plotted at distance, z — 1.34, such that the afterglow luminosity on day 
1 is comparable to ordinary GRBs (see also § 4). The SI observations are plotted with an 
assumed distance of 35.6 Mpc, the distance to the supernova. The subtraction of the known 
fluxes of sources Sib and Sic on day 1281 from the four BeppoSAX points reduces their 
values by ~ 4%, which is less than the size of the symbols used in the plot. We chose instead 
to plot the sum of all three SI sources in Figure 4 (Slgum) to indicate the "expected" flux 
from extrapolation of the BeppoSAX X-tslj light curve of SI and compare it with the flux of 
the Sla (the SN 1998bw) point only 

In contrast to the curve for SI or ordinary GRBs (Figure 3), the curve for S2 shows 
random variability and is nearly flat (assuming S(3') for the CXO flux of S2 on day 1281). 
It is consistent with a collection of variable X-ray sources, probably distant AGNs, whose 
sum sometimes exceeds and at other times falls below the BeppoSAX threshold. The total 
flux of the CXO observations of the sources within the larger (3' radius) error circle of S2 is 
within a factor of 1.5 of the brightest flux ever detected by BeppoSAX for S2 (day 1) and 
the two upper limits given by BeppoSAX on days 2 and 200. If the GRB occurred within 
S2 it either declined more rapidly than any other GRB ever studied before, in which case 
the observations of S2 offer no supporting evidence for the connection, or it created a most 
unusual afterglow that has not declined, in over three years. Moreover, the afterglow on day 
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1281 would have a luminosity orders of magnitude greater than other GRBs after day 50 
(Figure 3). The simplest conclusion is that S2 did not contain GRB 980425 and that the 
BeppoSAX detection was a collection of variable background sources. 

The hght curve of SI, on the other hand, shows a gentle decline to day 200, followed by 
a rapid fading, by factor of about 12 to 1.1 x 10^^ erg s~^, by day 1281 (here we compare 
the BeppoSAX \a\\ic at day 1 to Slgum)- This last data point is consistent with the light 
curves of other particularly luminous supernovae, e.g., SN 1993 J, that may have had high 
mass loss rates, but one must take care because SN 1998bw was a Type Ic supernova, which 
presumably occurred in a Wolf-Rayct star. Such stars are known to have a high wind velocity, 
and hence a low circumstcllar density. Bregman et al. (2003) have recently discussed the 
X-ray emission of young supernovae and find that most have an X-ray luminosity in the 
0.5 — 2 keV energy band less than 2 x 10^^ erg s~^. They do point out exceptional cases - 
SN 1978K, SN 1996J, SN 1998Z, SN 1995N, and SN 1998S - that have X-ray luminosities 
from 10^^-^ to 10^^ erg even a decade after the event, but these were all Type II. Given 
the exceptionally large kinetic energy inferred for SN 1998bw (Iwamoto et al. 1998; Woosley 
et al. 1999), perhaps it is not surprising that its luminosity after three years should place, 
e.g., about an order of magnitude above common Type Ic supernovae like SN 19941. 

The simplest conclusion here is that the brilliant emission of Sla during the first days 
was the X-ray afterglow of the relativistic ejccta that made GRB 980425. By day 1281 
however, we were seeing the energetic, but not especially relativistic ejccta of SN 1998bw 
colliding with the presupernova mass loss of its progenitor star. This hypothesis is discussed 
further in 5 5. 



5. Theoretical Interpretation 

The X-ray afterglow emission of ordinary GRBs is generally attributed to synchrotron 
emission from shocks as the blast encounters the interstellar or circumstcllar medium. Some 
useful scaling relations for blast waves in which each particle emits a fixed fraction, e, of the 
energy it gains in the shock, have been given by Cohen, Piran & Sari (1998). 

L{t) OC t-^^rn-2,)/{m+l)]-l 

with 

f2 + llr + 9 

3-e • ^2) 

Fundamentally, < e < 1 and 3 < m < 12, so that L oc to r^^/i^ ^ ^-1-69 xhis 
expression is for constant density. Ghevaher & Li (2000) also give expressions for the power 
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law scaling of afterglow light curves and find for a medium with p oar ^, 

Lit) oc r" (3) 

with a 1.75 to 2.17 for radiative blast waves and 1.38 to 1.75 for adiabatic blast waves if the 
index of the electron power distribution, p, is between 2.5 and 3.0. The curve L = lO^^t"^'^^ 
erg is plotted in Figure 3, and provides a reasonable description of the early X-ray 
afterglow lightcurve (when most of the energy is emitted at these frequencies). At observer 
times longer than a week, the blast wave would, however, be decelerated to a moderate 
Lorentz factor, irrespective of the initial value. The beaming and aberration effects are 
thereafter less extreme. If the outflow is beamed, a decline in the light curve is expected 
at the time when the inverse of the bulk Lorentz factor equals the opening angle of the 
outflow (Rhoads 1997). If the critical Lorentz factor is less than 3 or so (i.e. the opening 
angle exceeds 20°) such a transition might be masked by the transition from ultrarelativistic 
to mildly relativistic flow, so quite generically it would be difficult to limit the late-time 
afterglow opening angle in this way if it exceeds 20°. For reasonable conditions then, the 
power law declines of both XRFs and luminous GRBs, given in Figure 3, are what might 
be expected from very relativistic ejecta slowing in either a constant density medium or a 
circumstellar wind. 

However, the very slow decline of SN 1998bw needs a different explanation. GRB 980425, 
or at least that portion directed at us, was very weak. The total energy of its relativistic 
ejecta has been estimated as no more than 3 x 10^° erg (Li & Chevalier 1999) and most of 
that probably was not directed toward us. Wieringa, Kulkarni & Frail (1999) discuss the 
possibihty of two shocks associated with SN 1998bw, a relativistic one from the GRB and 
a slower moving, more powerful shock from the supernova, but until very late times when 
the blast reaches the termination of the presupernova wind, the thermal emission from the 
supernova itself is expected to be weak because of the low density of the wind (Chevalier 
2000). 

Two possibilities can be considered. First, that GRB 980425 was an "ordinary" (or 
somewhat subluminous) GRB observed off- axis. This has been suggested many times (e.g. 
Woosley, Eastman & Schmidt 1999; Nakamura 1999; Granot, Panaitescu, Kumar & Woosley 
2002; Yamazaki, Yonetoku & Nakamura 2003; Zhang, Woosley & Heger 2003) with different 
underlying assumptions regarding the angular distribution of the ejecta. Nakamura (1999) 
assumes that the jet has sharp edges and the peripheral emission comes from scattering. 
Woosley et al. (1999) assume that there is a distribution of ejecta energies and Lorentz 
factors and that, during the burst, we see only the low energy wing moving toward us. The 
other possibility is that GRB 980425 was deficient in energetic gamma-rays at all angles. This 
is not incompatible with the fact that it may have ejected 3 x 10^° erg of mildly relativistic 
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(r > 2) material, but concerns only the very relativistic ejecta, F > 200, thought responsible 
for harder GRBs. 

The X-ray light curve can help to distinguish these two possibilities. First, energetic 
as it may be, at early times the X-ray emission of the underlying subrelativistic supernova 
is probably neghgible. Though the mass loss rate of the Wolf-Rayet star progenitor may 
have been high, the wind velocity was also large implying a low circumstellar density and 
inefficient supernova emission (Chevalier 2000). Type lb and Ic supernovae are typically 
weak X-ray emitters (Figure 3; Bregman 2003). If we attribute the X-ray emission during 
the first 200 days to a weak, but on- axis GRB, then the X-ray afterglow should have faded 
with a power law not too different from the other GRBs. The fact that its decay is nearly 
flat is inconsistent with any relativistic blast wave in which the electrons emit a constant 
fraction of the energy gained in the shock (Cohen, Piran & Sari 1998), even in a constant 
density medium, and argues against an explosion with a single energy and Lorentz factor 
seen pole on. 

The data may be more consistent with a powerful burst seen off-axis (Granot et al. 
2002). Even along its axis, the burst may have been weaker than most, but the energy per 
solid angle could still have been orders of magnitude greater than along our line of sight, 
accounting for most of the ~ 3 x 10^° erg inferred from the radio (Li & Chevalier 1999). 
Because of relativistic beaming, initially we see only the low energy material moving toward 
us, but as the core of the jet decelerates, its afterglow is beamed to an increasing angle so 
that more and more energetic material becomes visible along our line of sight. Depending 
on the geometry, the afterglow luminosity could even temporarily increase. Granot et al. 
(2002) consider the appearance of GRB 980425 at various angles and conclude that the 
viewing angle needs to have been zSOg with Og the half angle of the most energetic part of 
the jet, otherwise the optical afterglow would have contaminated the supernova light curve 
unacceptably. One expects that the X-ray light curve would look similar to some of the 
plots of Granot et al. with the critical addition of low energy wings of ejecta as calculated 
by Zhang, Woosley & Heger (2003). This material would raise the luminosity at early times 
when almost nothing is seen of the central jet. A calculation of this sort must be a high 
priority for the community. 

As time passes, beaming becomes less important and the entire decelerating jet becomes 
visible, followed a little later by the underlying supernova. The light curve should then 
decline, as it did between days 200 and 1281 in Figure 3. One can estimate the time scale 
for when beaming becomes unimportant. For a circumstellar density distribution p = Ar^"^ 
Waxman, Kulkarni & Frail (1998), assuming typical mass loss and GRB parameters, estimate 
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a radius 



and a corresponding time 



tNR = ^ = 1.9(1 + .)^ yr (5) 



when the explosion has swept up a rest mass comparable to the initial relativistic ejecta. 
Here, A^, — A/b x 10^^ g cm~^, corresponding to a fiducial mass loss rate of 1.0 x 10"^ Mq 
y^-i ^Qj, ^ wind velocity of 1000 km/s, and E52 is the relativistic energy in 10^^ erg s"^^. 
One thus expects that for times of order several years the relativistic energy will have been 
radiated away and the emission will become isotropic. We can define this as the onset of the 
supernova stage. 

Evidence of a close association of the early X-ray emission to the overall GRB phe- 
nomenon is also presented by Bcrgcr, Kulkarni & Frail (2003). They have studied a sample 
of 41 GRB X-ray afterglows and they find a strong correlation between their X-ray isotropic 
luminosities (I>x,iso) (normalized to i = 10 hours after trigger) and their beaming fractions. 
We plot -Lx,iso as a function of the GRB isotropic equivalent 7— ray energy, E^so (Figure 5). 
The data are taken from Berger, Kulkarni & Frail (2003) and we have added data for GRBs 
031203, 030329 and 980425, as indicated on the plot. Several points are striking in this 
empirical relation: GRBs 031203 and 980425 fall well within the overall correlation and ex- 
tend the association by roughly three orders of magnitude. We have fitted the data without 
(dashed line) and with (sohd fine) all the outhers (980425, 031203,000210, 990705), with a 
power law of index 0.61 and 0.72, respectively; we conclude that the data are consistent with 
a trend extending roughly six orders of magnitude in X-ray luminosity. Parenthetically, the 
two outhers in the plot correspond to GRB 000210 (a 'dark' GRB) and GRB 990705 (a very 
bright GRB), and indicate that it may be possible to distinguish GRB subclasses by simply 
using their X and 7— ray properties, as also pointed out by Berger, Kulkarni & Frail (2003). 

Note that the convergence to a "supernova" at three to ten years here does not require 
all XRFs and GRBs to have a bright optical supernova following the GRB. Bright optical 
emission is a statement about the radioactivity the supernova made. X-ray emission at 3 
years is about its kinetic energy. If the variation in GRB energy and X-ray light curve is 
simply an effect of viewing angle, then one expects the empty parameter space in Figures 3 
and 5 to fill in with future observations. The recent GRB 031203 (Figures 3, 5) may well be 
the first of these 'gap' events (Rodriguez- Pascual et al. 2003; Soderberg et al. 2003a) if the 
reported redshift oi z = 0.105 (Prochaska et al. 2003) is confirmed. 
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6. Conclusions 

Our study resulted in the detection, on day 1281, of multiple X-ray point sources in the 
two BeppoSAX error boxes SI and S2 originally given (Plan et al. 2000) for the variable 
X-ray counterpart to GRB 980425 (§ 3.1). Based upon the accurately known radio location, 
the source Sla is definitely the supernova. The sum of all sources in S2 on day 1281 is similar 
to the maximum observed for this error box by BeppoSAX by a factor that is smaller than 
1.5. This is consistent with the hypothesis that S2 did not contain the GRB, but instead 
some variable X-ray sources that have not disappeared. SI, however, has been consistently 
observed to decline for ~ 1300 days. We conclude that the source SI contained the GRB 
and that SN 1998bw and GRB 980425 were the same explosion. 

Additional insight comes when we compare the X-ray hght curve of GRB 980425 with 

the broader family of high energy transients (§ 4). Comparison of published X-ray light 
curves for GRB, XRF and SNR supports a unification hypothesis^ similar to that proposed 
for AGN. Under this hypothesis all of these sources are associated with standard supernova 
explosions of massive stars. We can also distinguish a strong unification hypothesis in which 
the properties of the external medium are also standardized. There is considerable, though 
still inconclusive, evidence that the strong unification hypothesis is false, though it is still 
worth testing. These explosions are conjectured to produce an anisotropic, beamed com- 
ponent associated with a decelerating, ultrarelativistic outfiow and an unbeamed, isotropic 
component associated with the slowly expanding stellar debris. The flux associated with 
the beamed component depends upon the observer direction and declines rapidly with ob- 
server time; the closer to the axis the larger the flux and the more rapid is the onset to 
the typical afterglow decline. As the beamed component decelerates to become entirely 
non-relativistic, the observed flux will become independent of orientation. Under the strong 
unification hypothesis a one-parameter (inclination) family of X-ray light curves will be pro- 
duced, converging asymptotically to a single variation when the beamed component becomes 
non-relativistic. If only the unification hypothesis is true, we should still be able to observe 
these trends and estimate the inclination. Indeed, this is roughly what we observe when we 
plot the isotropic luminosities of GRB, XRF and SNR on a common scale (Fig. 3). It appears 
that after three years all explosions are subrelativistic, with X-ray luminosity dominated by 
the stellar debris, ~ 10^^ erg s~^. We therefore tentatively identify GRB, XRF and SNR as 
similar objects observed with small, medium and large inclination respectively. More specific 
to this paper, the observation that SN 1998bw and GRB 980425 follow a smooth light curve, 
which fits this pattern supports the claim that they are the same source. 

In § 5 we discussed two possible interpretations of these light curves based upon either 
a standard phenomenon viewed at different angles or explosions that eject variable amounts 
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of relativistic ejecta. We conclude that the observations, especially the slow initial decline 
rate, are more consistent with the "off-axis model" in which GRB 980425 was a much more 
powerful GRB seen at an angle greater than about three times the opening angle of the 
central jet. Emission at early times does not come from this central jet but from wings of less 
energetic material. After about three years the emission of all these high energy transients 
becomes isotropic and we see the sub-relativistic ejecta of the supernova interacting with 
the circumstellar wind. Thus all high energy transients have a common luminosity at three 
years due to their non-relativistic ejecta, but they follow different decay rates, depending 
upon the viewing angle, to reach there. 

Further, we discussed the stellar environment in that region of the galaxy ESO 184— G82 
where the supernova occurred. The supernova is one member of an X-ray doublet, both of 
which seem to be in the galaxy and to have spend part of their life cycles as ULXs. The 
projected distance between the two X-ray sources is ~ 300 pc, which is suggestive of some 
sort of a very active star forming region. 

Finally, we would like to stress that the existence of a relation between the dechne 
rate of the X-ray light curve during the first few weeks and its brightness implies that such 
measurements might be useful for diagnosing the character and subsequent evolution of a 
given high energy transient as well as constraining its distance. However, the sparse coverage 
of the current X-ray afterglow data does not allow us to address fundamental questions such 
as: Did the rapid decline of GRB 980425 continue? Did it/will it level out above or in 
the vicinity of other Type Ib/c supernovae? At what point does SN 1998bw become like 
ordinary supernovae? The current data enable us to make a prediction on the unification of 
the GRB-SN phenomena, which can only be vindicated with further observations, obviously 
of SN 1998bw, but also of the nearest XRFs and GRBs to fill in the missing parameter space. 
We strongly encourage, therefore, foUowup observations of nearby GRBs and XRFs for as 
long as the available instrumentation allows; we also encourage the calculation of off-axis 
models of the X-ray light curves, especially for a variable distribution of Lorentz factors and 
energies. 
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Fig. 1. — (a): Digitized Sky Survey (DSS) image of the area around GRB 980425. The outer 
(partial) circle represents the 8.0' error region for GRB 980425 from the BeppoSAX WFC. 
The full solid circles represent the 1.5' BeppoSAX NFI SI and S2 error regions. The dashed 
circle (radius= 15") within SI identifies the area covered by the GRB 980425 host galaxy 
ESQ 182 - G82. {b):CXO image of SI. The boxed region containing Slab (lower left) and 
Sic (upper right) is shown in (d). {c):CXO image of S2. Almost the entire error region (over 
97%) was covered by the CXO CCDs S2 and S3. Source S2d falls into the gap between these 
two CCDs, (d): CXO blow up of the square box in (b). The dashed circle (radius= 15") 
indicates the host galaxy area (also identified in (a)). A close-up view of the solid box shown 
here is seen in (e). {e):CXO blow up of the square box in (d). Sources Sla and Sib are 
clearly resolved. Sla is the source to the south and coincides with the established radio 
location of SN 1998bw. The two solid circles represent the la error regions for these sources 
with (r = 0.3") derived with comparative astrometry. {f):HST/STIS image of (e) taken in 
the 50CCD/Clear mode (no filter). The white circles are the CXO error circles from (e). 
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Fig. 2. — Light curves (left column; binned in 3000 s wide bins) and spectra (right column; 
binned in 0.232 kcV wide bins) of the five sources with more than 20 counts total within 
SI and S2 collected with CXO (0.3 — 10.0 keV). Notice source Sib is highly variable during 
the 50 ks observation. All other sources are either too faint to determine any variability or 
consistent with a constant persistent emission. 
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Fig. 3.— Compilation of GRB, XRF, SN I and SN II X-ray light curves (0.3 - 10.0 keV) 
presented as (isotropic) luminosity distances using a cosmology with Qm = 0.27, = 0.73, 
and Hq = 72 km s— 1 Mpc— 1. The XRF luminosities are calculated assuming two redshifts, 
z = 1 (solid lines) and z = 0.251 (dashed lines). The solid long line corresponds to a temporal 
decay of lO''^ tiavs^^'*' discussed in the text. 
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Fig. 4. — The X-ray light curves of sources SI and S2. The upper curve gives the two 
BeppoSAX observations and upper hmits for the error box of S2 during the first 200 days 
after the GRB trigger. The last points (day 1281) on that curve are the sum (S2sum) of all 
point sources observed with CXO and the sum (S2(3')) of all sources observed with CXO 
within an error circle of 3' radius. A distance of 2; = 1.343 has been assumed for all S2 points 
(see text). The arrow on day 1281 reflects the CXO detection hmit of a source placed at 
z — 1.343 or less. The lower curve is the light curve of source SI assuming a distance equal to 
that of ESO 184:— G82 (35.6 Mpc). The first four points are the BeppoSAX ohseivations; the 
fifth (upper) point at day 1281 is the sum of all SI CXO sources (Slsum)- The point below 
Slsum is the SN 1998bw luminosity after we subtract the fiux contribution from sources Sib 
and Sic from Slsum- 
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Fig. 5. — Isotropic equivalent luminosity of GRB X-ray afterglows, Lx, scaled to t = 10 
hours after the burst as a function of their isotropic equivalent 7— ray energy, Eiso. The data 
are taken from Berger, Kulkarni & Frail (2003); here we have added data for GRBs 031203, 
030329 and 980425, as indicated on the plot. The solid line is a fit including all data; the 
dashed hne is a fit exluding GRBs 980425, 031203, 000210, and 990705. 



